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ABSTRACT 

The period-luminosity sequences and the multiple periods of luminous red giant 
stars are examined using the OGLE III catalogue of long-period variables in the Large 
Magellanic Cloud. It is shown that the period ratios in individual multimode stars are 
systematically different from the ratios of the periods at a given luminosity of different 
period-luminosity sequences. This leads to the conclusion that the masses of stars at 
the same luminosity on the different period-luminosity sequences are different. An 
evolutionary scenario is used to show that the masses of stars on adjacent sequences 
differ by ^16-26% at a given luminosity, with the shorter period sequence being more 
massive. The mass is also shown to vary across each sequence by a similar percentage, 
with the mass increasing to shorter periods. On one sequence, sequence B, the mass 
distribution is shown to be bimodal. It is shown that the small amplitude variables on 
sequences A', A and B pulsate in radial and nonradial modes of angular degree i=0, 1 
and 2, with the i=l mode being the most common. The stars on sequences C' and C 
are predominantly radial pulsators (-^=0). Matching period ratios to pulsation models 
shows that the radial pulsation modes associated with sequences A', A, B, C' and C 
are the 4th, 3rd, 2nd and 1st overtones and the fundamental mode, respectively. 

Key words: stars: AGB and post-AGB - stars: oscillations - stars: variables: general. 


1 INTRODUCTION 

It is now well known that luminous red giants are variable 
and that the periods of the observed variability lie on dis¬ 
crete sequences in period-luminosity diagrams. [Wood et al.| 
119991 and Wood (20001 used data from MACHO observa¬ 


tions of stars in the Large Magellanic Cloud (LMC) to de¬ 
fine five sequences A, B, C, D and E in the (Iw, logP) and 
(K, log P) planes, respectively (where Iw is a reddening-free 
Weisenheit index defined as /-1.38(U-/)). The sequences A, 
B and C were identified as containing pulsating stars, while 
sequence E was found to consist of binary stars and sequence 
D was the sequence formed by the long secondary periods 
(LSPs) that occur in 30-50% of luminous pulsating red gi¬ 
ants. In this paper, the pulsating stars will be discussed. 
These stars lie on sequences A, B and C and two other se¬ 
quences A' and C' which are described below. Further ex¬ 
amination of the ellipsoidal variables on sequence E can be 


found in Soszyriski et al. (2004bh, Nicholls, Wood, & Cioni 


(20101, Nicholls & Wood (Wm and Nie & Wood (20141 


while studies of the origins of the LSPs can be found in Hin¬ 


kle et al. (20021, Olivier & Wood (20031, Wood, Olivier, & 
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Kawaler| (|2004[), [Nicholls et al.| ( |2009| ), jStother^ ( |2010| ) and 
Takayama, Wood, & Ita (20151. 

Following the discovery of sequences A, B and C, |Ita| 


et al. (20041 used extensive K band photometry combined 


with MACHO periods to show that sequence B was in fact 
two sequences which were renamed B and C' . The next 
refinement to the luminous sequences was by jSoszyiiski etj 
al. (2004a I who used OGLE data to show the existence of 
a new sequence on the short period side of sequence A. 


Soszyhski et al. (2004a 1 labelled this sequence a ,4 but in this 

paper the designation A' will be used (e.g. 

Derekas et al. 

2006||Tabur et al.|2010 

[Soszyhski, Wood, & Udalski^2013p. 

Soszyhski et al. 

(2004a 

1 also made the important discovery 

that sequence 1 

t consisted of three closely-spaced parallel 


sequences. jSoszyiiski et al.] ( |2007| ) claimed that sequences 
B and A' also consisted of three closely-spaced parallel se¬ 
quences (hereinafter referred to as sub-sequences). More re¬ 
cently, jSoszyhski fc Woodj ( |2013| ) have drawn attention to a 
probable faint sequence designated F between sequences C' 
and C. 


Wood et al. (1999 \ assumed that all the variable lu¬ 


minous red giants found by the MACHO observations were 


asymptotic giant branch (AGB) stars. However, Ita et al. 
(20021 and Kiss & Bedding (20031 showed that both AGB 


and red giant branch (RGB) stars were variable and that 
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2 P. R. Wood 


there was a slight offset in period between the sequences of 


RGB and AGB stars. Soszyhski et al. (2004a I also noted this 


offset but they did not find evidence for two slightly offset 
sequences for sequence A' which they suggested consisted 
of AGB stars only. However, this result was challenged by 


Tabur et al. (20101 who showed that an RGB component of 
sequence A' does exist. 

Various attempts have been made to identify the pul¬ 
sation modes associated with the pulsation sequences G, C', 
B, A and A' and the sub-sequences of B, A and A'. |Wood| 
et al. (1999 \ (see also Wood fc Sebo|[r996 \ showed that the 
radial fundamental mode could fit sequence C with reason¬ 
able assumptions about the mass of the star and that the 
first few radial overtone modes had periods similar to those 
of sequences B (which now includes G') and A, although 
the data available at the time was not extensive enough for 
a detailed matching of these latter sequences with modes. 
They also noted that the large-amplitude Mira variables fell 
on sequence C, making them fundamental mode pulsators. 
Additional evidence in favour of this modal assignment is 
provided by the fact that nonlinear fundamental mode radial 
pulsation models reproduce the light and velocity curves of 


local Mira variables quite well (e.g. Bessell, Scholz, & Wood 


1996 Ireland, Scholz, & Wood 20081. Recently, using ex¬ 


tensive observational data from OGLE, |Soszyhski fc Wood| 
(20131 showed that the period ratios in stars of sequences C 


and G' were consistent with those expected for fundamental 
mode pulsation of sequence G and first overtone pulsation 
on sequence C'. Takayama, Saio, & Ita (20131 compared 


period ratios for OGLE variables near the tip of the RGB 
with models and concluded that sequences G', B and A cor¬ 
responded to the radial first, second and third overtones. 


consistent with the assignments by Wood et al. (19991 and 


Soszyhski & Wood (20131. 


In contrast to the above results, Soszyhski et al. (20071 
concluded that the sequences A, B and C' corresponded 
to the radial 2nd overtone, first overtone and fundamental 
modes, respectively, when they used used radial pulsation 
models to try fit the sequences in the period-luminosity 


(PL) plane at RGB luminosities (see also Dziembowski & 


Soszyhski 20101. This means that the longer period sequence 


C containing the Miras has no explanation in terms of radial 


pulsation. This study was extended by Mosser et al. (20131 


who matched the A and B sequences at RGB luminosities 
to sequences of Kepler radial mode pulsators in a PL plane 
where frequency data was used as a proxy for luminosity. 
This result produces a close match between sequences B and 
A and the radial first and second overtones, respectively, in 


et al. (19991, Soszyhski & Wood (20131 and Takayama, Saio, 


agreement with Soszyhski et al. (20071 but not with Wood 


& Ita (20131. There thus exists two extant assignments of 


the radial modes to the PL sequences A, B, C' and C that 
differ by one radial order. 


In another study using Kepler data, Stello et al. (20141 
showed that in the luminous Kepler red giants, triplets of 
modes with periods corresponding to spherical degrees i=0, 
2 and 1 (in increasing order of period) are excited for a given 
radial order. The mode with l=\ has the largest amplitude 
for radial orders n=3 or more (the second or higher radial 
overtones). This suggests the OGLE sequences B, A and A' 
may be dominated by i=\ modes. 

In this paper, the period and luminosity data given in 
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Figure 1. The sequences of variable luminous red giants in the 
(WjKjlog PI) plane, where PI is the period of the mode with the 
largest amplitude and Wjk = K — 0.686(J — K) is a reddening- 
free measure of the luminosity. Mira variables are shown as red 
points (in the online version of the paper). The ellipsoidal vari¬ 
ables belonging to sequence E are not shown. 


the OGLE III catalogue of long-period variables (LPVs) in 
the LMG ( Soszyhski et al.|2009 1 is analysed in detail. Purely 
empirical results related to the modes, masses and evolution 
of the stars on the sequences are derived. Then the data are 
compared to pulsation models in order to throw further light 
on the modes of pulsation involved. 


2 THE SEQUENCES OF LUMINOUS RED 
GIANTS 

The sequences of luminous red giants in the LMG are shown 
in Figure The period data for this figure come from the 
OGLE III catalogue of LPVs in the LMG ( Soszynski et al. 
|2009p . All the variables in the catalogue are included in the 
analysis in this paper - no selection is made according to 
the catalogue types Mira, SRV or OSARG. In the cata¬ 
logue, three periods are given for each star with the first (or 
primary) period (PI) having the largest amplitude and the 
secondary periods P2 and P3 having progressively smaller 
amplitudes (most luminous red giants exhibit multiple peri¬ 
ods of variation). The luminosity in Figureis represented 
by WjK = K — 0.686{J — K) where the J and K magnitudes 
are from the 2MASS point source catalogue ( |Skrutskie et al.| 
20061. The individual sequences are labelled in Figure [^al 


though sequence E (produced by the binary ellipsoidal vari¬ 
ables) is omitted as these stars are not included in the cat¬ 
alogue. The tip of the RGB is clearly seen at Wjk ~ 11.3. 

The large-amplitude Mira variables are shown as red 
points in Figure]^ They clearly belong to sequence G. As 
noted in the Introduction, there are several lines of evidence 
which suggest that stars on sequence C are fundamental 
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Red giant pulsation modes, masses and evolution 3 


mode radial pulsators, although some authors have provided 
arguments against this modal assignment. In Section we 
will present what we believe to be compelling evidence that 
the stars on sequence C are indeed fundamental mode pul¬ 
sators. For these reasons, and to simplify some of the dis¬ 
cussion, in this paper it will be assumed from here on that 
the sequence C stars, including the Mira variables, are fun¬ 
damental mode radial pulsators. 

It is worth considering at this point what causes the 
existence of the sequences in diagrams like Figure Firstly, 
remember that the period used in Figure is the primary or 
largest-amplitude period. At a given luminosity, a star will 
have possible modes of pulsation with periods that could 
fall on any of the sequences or (as shall be shown below) 
between the sequences. The actual position of a star in the 
figure is determined by the period of the mode of largest 
amplitude. For self-excited modes, this mode is the one that 
has the strongest driving - see Soszyriski & Wood (20131 for 
a demonstration of this effect for 1st overtone pulsation on 
sequence C'. For modes that are excited stochastically by 
convective motions (solar-like oscillations), there is also a 
preferred pulsation period which is related to the acoustic 


cut-off frequency (e.g. Dziembowski et al. 2001 Mosser et 
^[M^lsTello et al.||2014|). 


The coloured points and short lines in the lower panel 
of Figure correspond to the ratios of the periods on the 
fit lines shown in the top panel i.e. they represent the ra¬ 
tios of the periods of the different period-luminosity se¬ 
quences. The horizontal position of the coloured point/lines 
is log P1[Wjk=12] for the line which is in the denomina¬ 
tor of the period ratio. The colour of the point/line is the 
colour of the line in the numerator. Coloured points corre¬ 
spond to period ratios of the fits to sequences A', A and B 
(these all have the same slope). For example, the red point at 
(logPl[WjK=12], logPi/PI) « (1.12, 0.17) represents the 
ratio of log PI at the middle of sequence A to log PI at the 
middle of sequence A', the red point near (1.3, 0.0) repre¬ 
sents this ratio for sequence A to itself and the red point near 
(1.52, -0.23) represents the ratio for sequence A to sequence 
B. Since the lines fitted to sequences C and C’ have a slightly 
different slope to that of sequence A, the period ratios of 
fits involving sequences C or C’ lead to the short coloured 
lines. For example, the short red line at (log P1[Wjk=12], 
logPi/Pl) ~ (1.65, -0.35) represents the ratio of log PI at 
the middle of sequence A to log PI at the middle of sequence 
C'. 

We now examine the structure in the lower panel of 
Figure in detail. 


3 PERIOD RATIOS 


Period-luminosity (PL) laws such as those in Figure [^pro¬ 
vide one type of diagram for the comparison of models and 
observations of variable stars. A more precise diagram for 
comparing observations and models is the Petersen diagram 
wherein period ratios in multimode pulsators are plotted 
against the longer of the two periods in the ratio. In this 
section, a modified Petersen diagram appropriate for LPVs 
is discussed. 

The upper panel of Figurej^shows the same PL diagram 
as Figure but now coloured straight lines have been fitted 
to each of the pulsation sequences C, C', B, A and A'. It has 
been assumed that the three sequences B, A and A' have the 
same slope dWjK/dlogPl = 4.444 while the two sequences 
C and C' have a slightly different slope dWjK/dlogPl = 
4.800. 


Since many of the variable LPVs are multimode pul¬ 
sators, they can be plotted in a Petersen diagram. The lower 
panel of Figure]^ shows a variant of this diagram where the 
logarithmic period ratios logP2/Pl and logP3/Pl in each 
star are plotted against the variable logPl[WjK = 12]. The 
latter quantity is the value of log PI obtained by project¬ 
ing the actual value of log PI along a line parallel to se¬ 
quence A until it meets the level Wjk = 12. Specifically, 
P1[Wjk = 12] = logPl -I- (Wjk — 12)/4.444. Stars with pri¬ 
mary periods on a given sequence are thus bounded in the 
same horizontal range in the lower panel of Figure with 
this range corresponding to the values of log PI for this se¬ 
quence when Wjk = 12. For LPVs, this form of Petersen 
diagram is much better at revealing preferred period ratios 
than the traditional Petersen diagram used by|Soszyhski et| 


al. (2004a I and Takayama, Saio, & Ita (20131 which plots 


Pshort/Pong against Pong. A form of Petersen diagram which 
used a quantity like P1[Wjk=12] on the horizontal axis was 
presented by Soszyhski et al.| (|2007[) . 


3.1 Observed period ratios for individual stars 
and for sequences 

The black points in the lower panel of Figure show ob¬ 
served period ratios for individual stars. The clustering of 
these points represents ratios of specific modes present in 
the data. Firstly, we start with the points in the vertical 
strip 1.2 < log P1[Wjk=12] < 1.4 which belong to stars 
with their primary period on sequence A. There is a large 
cluster of points near (logPl[WjK = 12], logPi/Pl) = (1.3, 
0.14) which one would expect to belong to stars on sequence 
A that show a secondary period on sequence B. Similarly, 
the cluster of points near (log P1[Wjk = 12], logPi/Pl) = 
(1.3, 0.3) should correspond to stars on sequence A that 
show a secondary period on sequence C' and a third cluster 
of points near (log P1[Wjk=12], logPi/Pl) = (1.3, -0.12) 
which should correspond to stars on sequence A that show a 
secondary period on sequence A'. However, the period ratio 
of sequences B, C' and A' to sequence A, represented by the 
green, magenta and orange points at log P1[Wjk=12]=1.3, 
respectively, do not coincide with these clouds of black 
points. It can be seen that this problem exists for stars with 
primary periods on any of the sequences. We will explore 
this situation in detail in Section |4] 

Another interesting feature of Figure is that there 
are clouds of points such as those near (log P1[Wjk=12], 
log Pi/Pl) = (1.45, -0.15) for which the primary period lies 
between the sequences. One can see in the upper panel of 
Figure that some stars do have their primary periods be¬ 
tween the sequences although they are not as common as 
the stars on the sequences. In fact, it turns out that the 
existence of stars between the sequences is related to the 
problem noted in the previous paragraph and we come back 
to this in Section [d] 


© 2014 RAS, MNRAS 000,[iHT6l 





































4 P. R. Wood 


8 


10 



12 


14 


1 


1.5 2 2.5 3 

log PI 




log(Pl[WJK=12]) 


Figure 2. Top panel: The black points show the stars in the (ITjki log PI) plane while the coloured lines are fits to the sequences. The 
coloured dots and thick coloured lines are projections of the coloured lines onto Wjk = 12 along lines of slope dVFjK/dlogPl = 4.444. 
Bottom panel: The black points show the logarithmic period ratios log P2/log PI and log P3/log PI in individual stars plotted against 
log P1[VFjk=12] (see text for the definition of this quantity). The coloured points and short lines correspond to the ratios of the periods 
on the fit lines shown in the upper panel (see text for full details). The colours can be seen in the online version of the paper. 


3.2 Closely spaced periods within a sequence 


A prominent feature in the lower panel of Figure is the set 
of sub-ridges first noted by Soszytiski et al. (2004a I and espe¬ 
cially prominent for stars on sequence A where five almost- 
horizontal sub-ridges can be seen in the cluster of points 


near (log P1[Wjk=12], log Pi/PI) = (1.3, 0.0). These five 
sub-ridges must be due to pulsation modes with very closely 
spaced periods. The spacing is sufficiently small that the in¬ 
dividual modes can not be distinguished in a PL diagram 
such as the top panel of Figure]^ The period separation of 
low-order radial modes in red giants is much larger than the 
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Table 1. Stars on the sub-sequences of sequence A 


Subgroup 

No. of stars 

Fraction of sequence A 

n(A) 

38526 

1.000 

n(0+l) 

639 

0.017 

n(0+2,2-1-1) 

1152 

0.030 

n(l+l,24-2,0+0) 

10866 

0.282 

n(l+2,2+0) 

2930 

0.076 

n(l+0) 

2754 

0.071 


period separation of the modes belonging to the five sub¬ 
ridges of sequence A so nonradial pulsation modes, as well 
as a radial pulsation mode, must be involved. 

The five sub-ridges of sequence A can be created 
by three closely-spaced modes each separated in log P by 
~0.022, provided that any of the three modes can be the 
primary mode in at least some stars. We assume here that 
this is the case (if one mode was always dominant, five modes 
would be required). There is some evidence that there are 
three sub-ridges, corresponding to two close modes, asso¬ 
ciated with sequence B (see the cluster of points around 
(logPl[WjK=12], logPi/Pl) = (1.53, 0.0)). This evidence 
will be made more apparent in Section 

As noted above, the three modes associated with se¬ 
quence A must include nonradial modes because of their 
close spacing. We assume that all stars on sequence A, and 
indeed all stars on any given sequence, have the same radial 
order of pulsation. Nonradial pulsation models of luminous 
red giants (e.g. Dziembowski et al.|2001 Mosser et al.|20f3 


Stello et al.|2014 l show that for modes of a given radial or¬ 

der but different angular degree £ = 0, 1 or 2 (the only ones 
likely to be observed in light curves), the £ = 1 mode has 
the longest period, the £ = 0 mode has the shortest period 
and the £ = 2 mode has an intermediate period. 

Figure provides a method for determining which of 
the three modes within sequence A is dominant. Among the 
stars with their primary period on sequence A' (the points 
with 1.05 < log P1[Wjk=12] < 1.2), those that also have 
the closely-spaced modes associated with sequence A can 
be seen near (log P1[Wjk=12], logPi/Pl) — (1.12, 0.12). 
There are three horizontal bands here with the same spacing 
as the sub-ridges associated with sequence A. The upper 
band is by far the most prominent. Assuming that one mode 
dominates sequence A' (because there are no prominent sub¬ 
ridges near (log P1[Wjk=12], logPi/Pl) = (1.12, 0.0)), this 
band structure suggests that it is the longest period mode 
on sequence A that generally has the largest amplitude. As 
noted above, this is the £ = 1 mode. The dominance of 
the £ — 1 mode found on sequence A in the present high 
luminosity red giants is consistent with the findings of|Stello| 


et al. (20141 from Kepler observations but different from 


that of Mosser et al. (20131 who find the radial £ = 0 mode 
becomes dominant at high luminosities in their analysis of 
Kepler stars. 

Estimates of the relative amplitudes of the three modes 
that make up the three sub-sequences of sequence A can be 
found by examining the populations of the five sub-ridges of 
sequence A in the Petersen diagram. An exploded view of 
these sub-ridges is shown in the top panel of Figure]^ Boxes 
have been drawn to enclose the sub-ridges and the sub-ridges 
have been named after the possible modes occurring in them. 


CO 

oi 

II 


Si 

o 
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log(Pl[WJK=12]) 



Figure 3. Top panel: The small section of the Petersen dia¬ 
gram containing the five sub-ridges formed by the three pulsation 
modes forming sequence A. The five boxes enclose the five sub¬ 
ridges and their names are shown to the left of the boxes. Bottom 
panel: The amplitude distribution of the £ = 0 and £ = 1 modes 
of sequence A when these two modes are present together. See 
text for further details. 


Under our assumption that the three modes involved are the 
£ = 0, 1 and 2 modes of a single radial order and that the 
period increases from £ = 0 to 2 to 1, a star found in the top 
box must have £ = 0 for its primary period PI and £ = 1 
for the secondary period (P2 or P3). This box has been la¬ 
belled 0+1. Similarly, the lowest box 1+0 must have £ = 1 
for its primary period PI and £ = 0 for its secondary period. 
The central box is made up of stars where the same mode is 
detected twice with slightly different periods and the mode 
involved could be any of the £ = 0, 1 and 2 modes. The 
two intermediate boxes are labelled by the possible pairs of 
modes that could be present. Note that the period separa¬ 
tion A log P ~ 0.043 of the £ = 0 and £ = 1 modes obtained 
from the positions of the sub-ridges agrees with the period 
separation of these modes given by the observations and the 
models of Stello et al. (2014l, giving added support for the 
mode identifications. 

The fact that the ridges in boxes 0+1 and 1+0 are not 
strictly parallel to the horizontal axis means that A log P 
varies systematically by a small amount across sequence A. 
Quantitatively, for stars on sequence A, A log P decreases by 
~0.0045 as log P1[Wjk=12] increases by 0.1. An explanation 
for this behaviour is given in Section [7.2| 

The number of stars in each box is given in Table 
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6 P. R. Wood 


The total number n(A) of stars with their primary period 
PI on sequence A is also given. This latter number is the 
total number of stars with log P1 [Wjk= 12] lying between 
the values corresponding to the two vertical lines in Figure]^ 
Stars that have both the i = 1 and £ = 0 modes detected in 
their light curves lie in boxes 1+0 and 0+1. It is clear that 
n(0 + 1) is much smaller than n(l + 0) which shows that 
the 1=1 mode is much more likely to be detected as the 
primary pulsation mode than the £ = 0 mode and this in 
turn suggests that the £ = 1 mode is likely to have a larger 
amplitude than the £ = 0 mode. This is shown directly in 
the bottom panel of Figure where histograms of modal 
amplitudes are shown (the amplitudes come from the OGLE 
III catalogue). The top sub-panel shows the amplitudes of 
the £ = 0 and £ = 1 modes for stars in box 1+0, the middle 
sub-panel shows the amplitudes for stars in box 0+1 and the 
bottom sub-panel shows the sum of the top and middle sub¬ 
panels. It is clear that the £ = 1 mode tends to have the 
larger amplitude in stars where the £ = 0 and £ = 1 modes 
are both present although the non-zero population in box 
0+1 shows that the £ = 0 mode can sometimes have a larger 
amplitude than the £ = 1 mode. 

The amplitude distributions shown in Figurej^apply to 
stars with two modes on sequence A. However, about 80% of 
sequence A stars have only one detected mode on sequence 
A (see the fractions in the last column of Table and note 
that stars on the middle sub-ridge have two detections of 
the same mode). Assuming that the relative amplitude dis¬ 
tributions shown in Figure]^ are typical, the stars with only 
a single mode on sequence A will be mostly £ = 1 pulsators, 
in agreement with the result derived above from the stars 
that have their primary mode on sequence A' and a sec¬ 
ondary mode on sequence A. 

One factor that has not been mentioned above is the 
orientation of the pole of the pulsating star relative to the 
observer in the case of nonradial pulsation. Different orien¬ 
tations will cause different apparent amplitudes of pulsation 
for a nonradial pulsator, with a zero apparent amplitude be¬ 
ing possible for specihc values of £, m and the orientation. 
As a result, the true numbers of ^ = 1 and £ = 2 nonradial 
pulsators relative to the number of ^ = 0 radial pulsators 
should be higher than indicated by the numbers in Table 


4 MULTIPLE MODES IN THE PL DIAGRAM 

In Section [3.1| it was shown that period ratios in individual 
stars do not agree with the period ratios of the mean se¬ 
quences. A vivid demonstration of this situation is shown in 
Figure Looking at panel A, it can be seen that the sec¬ 
ondary periods form three strips parallel to sequence A, but 
they are closer to sequence A than are sequences A', B or 
CL This situation, wherein the secondary modes are closer 
to the primary mode than are the adjacent sequences to the 
central sequence, exists for sequences A', A and C. It also 
exists for the modes shortward of the primary mode of se¬ 
quence B and for the mode longward of the primary mode for 
sequence C'. The only exception to this rule occurs for the 
mode longward of the primary mode of sequence B (which 
is further from sequence B than is sequence C') and for the 
mode shortward of the primary mode of sequence C' (which 


is further from sequence C' than is sequence B). We come 
back to this exception later. 

We now seek to explain why the green strips in Figure]^ 
do not align with the sequences. Firstly, it is assumed here 
that the pulsation modes that produce the green strips are 
the same pulsation modes as those that produce the cor¬ 
responding, neighbouring sequences (where the correspon¬ 
dence is obtained from the ordering in log P of the strips and 


sequences). This assumption will be justified in Section 5.1 


In this situation, the period of the same mode must be dif¬ 
ferent on a green strip and on its corresponding sequence. 
Now, in general, the period of a pulsation mode can be writ¬ 
ten as P oc /M^ where ot and /3 are both positive and 
for red giants typically range from 1.5-2.3 and 0.5-1.9, re¬ 
spectively (e.g. |Fox fc Wood 19821. In addition, the posi¬ 
tion of the giant branch at a given luminosity means that 
Tefi = Te{f{L, M). Combining these two relations with the 
definition L = yields P = P{L,M). Thus if a 

green strip and its corresponding sequence have the same 
mode but a different P at a given L then the masses of the 
stars in the green strip and its corresponding sequence must 
be different at a given luminosity. 

Finally in this section, we return to the exceptional 
cases of (1) the green strip originating from sequence B 
which lies on the long period side of sequence C' and (2) 
the green strip originating from sequence C' which lies on 
the short period side of sequence B. Remember that in Sec¬ 
tion it was deduced that the £ = 1 mode was the dom¬ 
inant mode in sequence A and that this mode has a period 
which is ~0.043 longer in logP than the radial mode with 
£ = 0. Recall also that there appear to be at least two modes 
(of the same radial order) associated with sequence B. We 
will show later that these are most likely the £ = 1 and 
£ = 0 modes, with the former being most common. Finally, 
as noted in Section]^ stars on sequence C appear to be pul¬ 
sating in the radial fundamental mode. With these facts in 
mind, the exceptional cases noted here can be explained nat¬ 
urally if stars on sequences A', A and B are most commonly 
pulsating in the £ = 1 mode in both their primary pulsation 
mode and their secondary modes P2 and P3, while stars 
on sequences C' and C are most commonly pulsating in the 
radial £ = 0 mode in both their primary pulsation mode and 
their secondary modes. Thus the stars with primary modes 
on sequence B will have logP for their secondary modes 
increased relative to logP of the radial mode pulsators on 
sequence C' by A logP ~ 0.043, thereby causing the ex¬ 
ceptions noted in points (1) and (2) above. Note that the 
offset A logP ~ 0.043 is derived from the populous RGB. 
Since the slope of sequences C and G' is steeper than the 
slope of sequences A', A and B, this offset will increase with 
luminosity and will be larger for AGB stars. 

An alternative explanation for the exceptional cases 
would be for the stellar mass to decrease from sequence A' 
to A to B, then to increase to sequence C' then decrease 
again to C. This seems contrived and is not consistent with 
the evolutionary picture described in the next section. 


4.1 The stars on sequence D 

Although generally not part of this study, in the last panel 
of Figurej^the secondary periods of stars with their primary 
period on sequence D are shown as green points. This plot 
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Figure 4. Plots of VPjk against the primary period PI for stars in the OGLE III catalogue (black points). Each panel is labelled by a 
sequence name and the stars with primary periods lying in a channel running along the centre of the named sequence have been plotted 
as orange points. Secondary periods P2 and P3 for the stars with PI in the orange channel are plotted as green points. For clarity, 
green points close to the orange channel are not plotted: the excluded points have | logP2 — log Pc | < 0.07 and | logP3 — log Pc | < 
0.07, where log Pc is logP at the centre of the channel. The excluded points are second detections of the mode represented by the orange 
points. The colours can be seen in the online version of the paper. 
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Figure 5. Plots of Wjk against the primary period PI for stars 
in the OGLE III catalogue (black points). Evolutionary channels 
cross each sequence from left to right and stars whose primary 
periods fall in these channels are plotted as orange points. The 
secondary periods P2 and P3 of the stars with PI in the orange 
channels are plotted as green points. See text for further details. 
The colours can be seen in the online version of the paper. 


highlights the fact that the periods of pulsation associated 
with the LSPs of sequence D lie mainly on or near sequence 

B. 

4.2 The stars on sequence F 

The faint sequence F between sequences C' and C is marked 
in Figure In panel C' of Figure it can be seen that 
secondary periods of stars whose primary period lies on se¬ 
quence C' form a sequence close to sequence F. This suggests 
that the sequence F stars are possibly fundamental mode 
pulsators like the stars on sequence C. However, it is not 
clear why the sequence F stars should form a distinct se¬ 
quence. One possible explanation is that they are high mass 
LPVs (so that their periods are shorter than the periods 
typical of sequence C stars) formed in a relatively recent 
burst of star formation. 


5 THE EVOLUTION OF LPVS 


5.1 Tracks in the PL diagram 

For most of the stars in PL diagrams such as Figure 
mass loss is small so that the evolutionary track of a star 
is represented by a line of constant mass and pulsation 
mode. An examination of pulsation models shows that in 
the (Wjk, logP) plane, lines of constant mass and pulsation 


mode have a smaller slope than the sequences (e.g. Soszyhski 
fc Wood|2013 1. In Figure]^ channels similar to those in Fig 
ure 1^ have been drawn but now the channels are roughly 
parallel to lines of constant mass and pulsation mode so 


they represent evolutionary tracks. The slopes of the tracks 
are estimated from the pulsation models described in Sec¬ 
tion For convenience of plotting, at lower luminosities 
on sequences A', A and B, the orange channels have been 
drawn more tilted relative to the sequences than the models 
indicate but this does not qualitatively affect the empirical 
evolutionary scheme presented here. 

Consider stars in Figure whose primary period lies 
in the orange strip through sequence A. As stars in this 
channel evolve to higher luminosity, they slowly move from 
the short period side of sequence A to the long period side. 
At the same time, there is a declining prominence of the 
secondary period corresponding to the next highest radial 
order (the green strip to the left of the orange channel in 
sequence A) and an increasing prominence of the secondary 
period corresponding to the next lowest radial order (the 
green strip to the right of the orange channel in sequence 
A). This means that the primary period most likely detected 
by observations will move to a lower radial order as the lu¬ 
minosity increases. Thus the orange points representing the 
evolving stars in Figure will jump from the long period 
side of sequence A to the short period side of sequence B 
because the highest amplitude mode has moved from that 
associated with sequence A to that associated with sequence 
B. Evolution will then continue up and across sequence B 
before a jump to sequence C' and so on. Note that in the 
empirical scheme shown here, the point of transition of the 
orange channels from sequence to sequence is judged by eye 
since current pulsation models do not give reliable theoreti¬ 
cal estimates of mode amplitudes. 

When a star reaches the RGB tip, it will transit to 
a lower-luminosity core-helium burning phase before rising 
again up the AGB with a slightly shorter period than it had 
at the same luminosity on the RGB. This situation is shown 
in Figure for the stars in the channel on sequence B. The 
AGB evolution is shown only after the RGB tip luminosity 
(Wjk ~ 11.3) is exceeded. 

5.2 Mass variations through the sequences 

The evolutionary scenario described in the previous sub¬ 
section has implications for the masses of stars on the se¬ 
quences. Note firstly that, on a given sequence, stars of dif¬ 
ferent mass will follow essentially parallel tracks but the 
higher mass stars will be found at higher luminosities. This 
means that the mass increases up a sequence, it increases 
to shorter periods across a sequence at a given luminosity, 
and it also increases from sequence to sequence at a given 
luminosity as the radial order increases (the sequence period 
decreases). Note that a mass difference between sequences 
was invoked above in order to explain why the period ra¬ 
tios between the sequences did not match the period ratios 
in individual stars. The evolutionary scenario provides an 
understanding for this empirical finding. 

The periods of stars in the various sequences formed 
in PL diagrams by using the primary or secondary periods 
can be used to estimate the mass variation from sequence 
to sequence and the mass variation across sequences. Fig¬ 
ure shows a histogram of the number of stars in a hor¬ 
izontal band defined by Wjk = 12 ± 0.25 in panel A of 
Figure]^ Here we are interested in exploring the mass vari¬ 
ation from sequence B to sequence A and the mass varia- 
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Figure 6. Histogram of the number of stars in a horizontal band 
defined by Wjk = 12 it 0.25 in panel A of Figure]^ The hori¬ 
zontal position of each star is given by the value of its coordinate 
logP[lFjK=12]. The black curve is obtained by using the pri¬ 
mary period PI of each star. The orange segment of this curve 
represents the sequence A stars in the orange channel through 
the middle of sequence A. The green curve is obtained by using 
the secondary periods P2 and P3 of those stars in the orange 
channel (i.e. this is the histogram of the green points in panel A 
of Figure]^. The colours can be seen in the online version of the 
paper. 

tion across sequence A. The positions of sequences A and 
B are marked on Figure The peak in the green distri¬ 
bution at logP[WjK=12] ~ 1.44 is formed by sequence A 
stars whose secondary periods are produced by the pul¬ 
sation mode which produces sequence B. The difference 
AlogP ~ 0.081 in the positions of this peak and the se¬ 
quence B peak is caused by the different masses of stars on 
sequences A and B. Pulsation models for low radial order 
modes at the luminosities and masses of the stars involved 
here show that, at a given luminosity, P oc ® (this 
result comes from the pulsation models described in Sec¬ 
tion]^. Thus the shift AlogP ~ 0.081 corresponds to the 
average mass of stars on sequence A being approximately 
26% larger than the average mass of stars on sequence B. 

Mass variations within a sequence are also required. 
These mass variations make a large contribution to the fi¬ 
nite width in logP of the sequences, with the mass increas¬ 
ing to shorter periods across the sequences. Let us consider 
sequence A. The total width of this sequence (represented 
by the full-width at half maximum) is AlogP ~ 0.098. 
As well as the mass variation across the sequence, there 
is also a contribution to the sequence width from the pe¬ 
riod variability represented by the vertical width of the sub¬ 
ridges in the top panel of Figurej^ This contributes approx¬ 
imately AlogP ~ 0.01 to the observed width of sequence 
A, leaving the contribution from mass variation alone at 
AlogP ~ 0.088. The sub-sequences of sequence A (i.e. the 
multiple modes) will also make a small contribution to the 
width but we ignore this. The width AlogP ~ 0.088 means 
that the mass increases by about 29% from the longer period 
to the shorter period side of sequence A. 

In the evolutionary scenario described above, the mass 
variation across a sequence should be equal to the mass 
variation between sequences. Remember that stars jump be¬ 


tween sequences from the long period edge of one to the short 
period edge of the other. Thus the masses of stars on the long 
period side of a sequence should be the same as the masses 
of stars on the short period edge of the adjacent longer pe¬ 
riod sequence. Consequently, the mass difference between 
the short period edges of two adjacent sequences will be the 
mass difference across the shorter period sequence. The re¬ 
sults given above for sequences A and B show that the mass 
difference between these two sequences is indeed close to the 
mass difference across sequence A. 

The procedure applied above to determine the differ¬ 
ence in the average masses of stars on sequences A and B 
can be applied to other pairs of sequences in which stars 
pulsate mostly in the same mode. We find that stars on se¬ 
quence A' are ~17% more massive than stars on sequence A 
at luminosities around Wjk = 12. The difference in the aver¬ 
age masses of stars on sequences C' and C cannot be readily 
determined at this luminosity since the mode of sequence C 
is rarely seen there. However, near Wjk = 10.5, which is on 
the AGB, many stars exhibit primary and secondary oscilla¬ 
tions in both of the modes associated with sequences C' and 
C. The period shifts indicate that the stars on sequence C' 
are ~16% more massive than those on sequence C. To de¬ 
rive this result we used the formula P oc M~^'^ rather than 
P oc which is appropriate at lower luminosities on the 

RGB (these formulae are derived from the pulsation mod¬ 
els described in Section]^. Finally, we consider sequences 
B and C'. Gomparing the masses of stars on these two se¬ 
quences is difficult at any luminosity because of the paucity 
of stars that exhibit both the sequence B and C' modes, and 
because the stars with primary periods on sequence B ap¬ 
pear to be £ = 1 pulsators while stars with primary periods 
on sequence C' appear to be ^ = 0 pulsators. We do not 
make any estimate of the mass difference between these two 
sequences. 


6 THE LUMINOSITY DEPENDENCE OF 

PERIOD RATIOS 

Although Figure [^demonstrates the existence of secondary 
periods in many LPVs, it does not indicate how the presence 
of secondary periods varies with luminosity. The luminosity 
dependence of the secondary periods can be seen in Figure 
For each sequence, this figure plots Wjk against the sec¬ 
ondary period relative to the primary period whose position 
is marked by a vertical blue line. The vertical strips of black 
points show where individual secondary modes are present 
in the stars. The radial order of the modes associated with 
these secondary modes are indicated by coloured numbers, 
with radial order 1 corresponding to the fundamental mode 
(see Sectionj^for the justification of the assignment of radial 
orders). The positions of the individual sequences are indi¬ 
cated by coloured dots in each panel. The fact that the dots 
and numbers of the same colour do not align horizontally 
(except for the reference mode) is yet another demonstra¬ 
tion of the fact that the ratio of the periods in individual 
stars do not match the ratio of the periods of the sequences. 

Panel A' of Figurej^shows the ratio of secondary to pri¬ 
mary period for stars whose primary period lies on sequence 
A'. Looking at this panel, we see three closely-spaced verti¬ 
cal bands under and slightly left of the red 4. These bands 
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Figure 7. H^jk plotted against the ratio of secondary period to the primary period for stars whose primary periods lie in a strip of 
width AlogPl = 0.12 down the centre of each sequence. Here logPj-ei Is defined as logPrei = log(P2/Pl) + logPcentre or log-Prel = 
log(P3/Pl) + logPcentre, where logPcentre Is logP at the Centre of the sequence when Wjk = 12. The blue vertical line in the panel for 
each sequence is situated at logPcentre- The coloured dots indicate logPcentre for each sequence using the colour coding for the sequences 
adopted in Figure]^ The coloured numbers count the observed radial orders in each panel and these numbers are aligned consistently 
from panel to panel, starting at 1 for the mode with the longest observed period in panel C. The colours can be seen in the online version 
of the paper. 
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are caused by the three modes occurring on sequence A. 
The longest period mode (the 1 = 1 mode) becomes domi¬ 
nant as the luminosity increases. This result is also seen for 
stars with their primary period on sequence A (see panel A). 
There are five vertical bands associated with sequence A and 
these are caused by the the three modes present within this 
sequence. As discussed in Section |3.2[ the shortest period 
band is produced by stars whose primary mode has 1=1 
and whose secondary mode has i = Q while the longest pe¬ 
riod band is produced by stars whose primary mode has 
I = Q and whose secondary mode has i = 1. The shortest 
period band extends to higher luminosities than the longest 
period band which suggests that the 1 = 1 mode is more 
prominent as the primary period (i.e. the one with largest 
amplitude) at higher luminosities. Once again, we note that 
the dominance of the 1=1 mode at higher luminosities on 
sequence A is consistent with the findings of |Stello et al.| 
(20141 from Kepler observations but different from the find¬ 


ings of Mosser et al. (20131 who suggest that the radial £ = 0 


mode becomes dominant at high luminosities. 

The stars with their primary mode on sequence B show 
thre^ vertical bands associated with the primary mode 
(panel B of Figure [^. This structure indicates the pres¬ 
ence of two close modes but they seem to be confined to 
RGB stars only (Wjk > 11.3). Only one mode appears to 
be present on the AGB. 

Panel B in Figure reveals the presence of apparent 
extra modes not present in any other panel. In a simple se¬ 
quential assignment of radial orders, there appear to be sec¬ 
ondary modes present corresponding to radial orders 5 and 
6. However, looking at panel A' immediately above panel B, 
it can be seen that the mode structure associated with the 
radial orders 4 and 5 in panel A' is reproduced exactly in 
panel B. Another point to notice is that the horizontal posi¬ 
tion of radial order 2 in panel A' coincides exactly with the 
position of sequence B. We conclude from these facts that 
stars whose primary mode lies on sequence B are of two 
types: (1) stars pulsating in the third radial order with sec¬ 
ondary periods in the second and fourth radial orders, and 
(2) stars pulsating in the second radial order with secondary 
periods in the third, fourth and fifth radial orders. The type 
2 stars are more massive than type 1 stars and they have 
masses typical of stars on sequence A'. These type 2 stars 
are really stars belonging to sequence A' but the observed 
amplitude of the radial order 2 pulsations is larger than the 
observed amplitude of the radial order 5 pulsations so they 
end up on sequence B. The apparent extra radial orders in 
panel B have been labelled by their true radial orders using 
numbers in brackets. 

If we consider the stars on sequence B whose primary 
mode is of radial order 3 (panel B of Figure]^, we see that 
the secondary modes pulsating in radial order 4 form two 
bands with the shorter period band being overwhelmingly 
dominant (the apparent band on the short period side of the 
radial order 4 group belongs to radial order 3 of the type 2 
stars examined in the last paragraph). Given that there are 


^ The apparent splitting near Wjk = 12 of the outer bands of 
sequence B are cause by the 1-year alias of the primary period, so 
this feature is not real. A similar 1-year alias can be seen in the 
bands of sequence A near W]k = 11-2 if the figure is expanded. 


two possible primary modes of radial order 3, the dominance 
of the shorter period band in panel B suggests that the dom¬ 
inant primary mode is the longer period one. The separation 
A log P « 0.04 of the two modes associated with radial order 
3 (derived from the separation of the vertical bands associ¬ 
ated with this order in panel B of Figure is consistent 
with the separation of the t = Q and the t = 1 modes of 
radial order 3 in the models of Stello et al.| ( |2014 1. These 
results therefore suggest that the t = Q and the 1=1 modes 
are present on sequence B and that the 1=1 mode is dom¬ 
inant. The direct finding that the 1=1 mode is dominant 
on both sequences B and A is consistent with the sugges¬ 
tion in Section 13 that the £ = 1 mode is dominant on each 
of sequences B, A and A'. 


7 COMPARISON OF OBSERVATIONS AND 
RADIAL PULSATION MODELS 

In this section, period ratios observed in individual stars are 
compared to period ratios in radial pulsation models of red 
giants. The linear, non-adiabatic pulsation models were com¬ 


puted with the code described in Wood & Olivier (20141 us¬ 


ing the same input parameters (abundances, mixing-length, 
turbulent viscosity parameter, core mass) as in that paper. 
Evidence was presented in Section that, in a given star, 
both the dominant primary pulsations and the secondary 
pulsations correspond to modes of the same spherical de¬ 
gree £. Thus, although the nonradial £ = 1 mode appears 
to dominate the primary pulsations in sequences A', A and 
B, provided the ratio of the ^ = 1 to the £ = 0 periods is 
approximately the same for radial orders 3 to 5 associated 
with these sequences, the observed period ratios in individ¬ 
ual stars should be the same as the period ratios in radial 
pulsation models. In Figures and the observed period 
ratios in individual stars are compared to the period ratios 
in radial pulsation models. 

In Figure]^ it is assumed that the radial mode associ¬ 
ated with sequence G (the Mira sequence) is the fundamental 
mode (radial order 1) and that sequences C', B, A and A' 
correspond to successively higher radial orders of pulsation. 
It can be seen that for each of the panels in Figure the 
model period ratios match the observed period ratios well for 
some mass in the range 0.8 to 2.4 M© except for the highest 
mass and luminosity models in panels G' and G. This prob¬ 
lem is discussed in Section [7.3| Note that in panel B, there 
is a mismatch between the models and those vertical bands 
that really belong to sequence A' (as described in Section]^ 
but this is expected because of the different masses of stars 
on sequences A' and B. Overall, these results provide strong 
evidence that the primary pulsation modes associated with 
sequences C, G', B, A and A' belong to radial orders 1, 2, 3, 
4 and 5. It also provides strong evidence that both primary 
and secondary pulsations in a given star tend to belong to 
the same £ value. 

In Figure it is assumed that the radial mode asso¬ 
ciated with sequence C' is the fundamental mode (radial 
order 1) and that sequences B, A and A' correspond to suc¬ 
cessively higher radial orders. Glearly, the mode marked by 
the cyan number 1 (sequence C) has no explanation in terms 
of radial modes with this assignment of radial orders. Fur¬ 
thermore, even for higher radial orders, the ratios of mode 
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Figure 8. Same as Figure[^with the addition of radial pulsation models showing the logarithmic period ratios for the first five modes. 
It is assumed that sequence C corresponds to the fundamental mode (radial order 1) and that sequences to A' correspond to the 
first to fourth overtone (radial orders 2 to 5), respectively. In each panel, the model period ratios are taken relative to the radial order 
appropriate for the sequence assigned to the panel. Models with masses of 0.8, 1.6 and 2.4 Mq are represented by red, orange and green 
lines, respectively. Both AGB and RGB models are plotted for M equal to 0.8 and 1.6 Mq (there is no RGB for M = 2.4Mq). The 
AGB and RGB models are only separable for the fundamental mode. The colours can be seen in the online version of the paper. 
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Figure 9. Same as Figure]^ but now it is assumed that sequence corresponds to the fundamental mode (radial order 1) and that 
sequences B, A and A' correspond to the first to third overtone (radial orders 2 to 4), respectively. Only the first four radial modes of 
the models are shown. 
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periods in the models fail to match the observed ratios. We 
therefore reject this modal assignment. 


7.1 Comparison with previous modal assignments 

The finding here that sequences C, C', B, A and A' corre¬ 
spond to radial orders 1, 2, 3, 4 and 5 is consistent with most 
earlier findings (e.g. [Wood &: Sebo||1996 Wood et al.||1999 


[Takayama, Saio, &: Ita||2013p . However, this assignment is 
inconsistent with the suggestion of Soszyhski et al. (20071 
(20131 that sequence C' corresponds to 


Mosser et al. 


and 

radial order 1. 

In terms of nonradial modes, we find direct evidence 
for the presence of the 1=1 mode on sequences A and B 
and for the 1 = 2 mode on sequence A. In their analysis of 
OGLE III data, Takayama, Saio, fc Ita| | |2013[ ) claim to find 
the presence of the two nonradial modes. Firstly, they asso¬ 
ciate the ^ = 2 mode of radial order thre^with sequence B 
while the only nonradial mode we find has 1 = 1. Secondly, 
they find the 1 = 1 mode of radial order 5 associated with se¬ 
quence A whereas we associate radial order 5 with sequence 
A' and we find the 1=1 mode of radial order 4 associated 


with sequence A. Because Takayama, Saio, & Ita (20131 only 


looked at LPVs on the RGB designated as OSARGS in the 
OGLE III catalogue, all LPVs with periods on sequences A' 
and C were excluded (by definition, OSARGS on the RGB 
do not have these modes) so they had an incomplete picture 
of the modes involved. In addition, the Petersen diagram 
they used meant that a number of modes overlapped on the 
figure causing confusion. 


7.2 The mass dependence of the period ratios 

In the lower panel of Figure the various clusters of points 
representing period ratios of specific pairs of modes tend 
to form ridges that are not horizontal. For example, we now 
know that the ridge at (log P1[ILjk=12], logPi/PI) « (1.3, 
0.3) corresponds to stars whose primary pulsation mode be¬ 
longs to radial order 4 and whose secondary pulsation mode 
belongs to radial order 2. This group of points, and most of 
the other groups of points, form ridges such that the period 
ratio becomes smaller when moving from the long period 
side of the primary mode sequence to the short period side. 
It was shown in Section [5.21 that the mass of stars in a se¬ 
quence increases when moving from the long period side to 
the short period side of the sequence. This suggests that the 
cause of the non-zero slope of the ridges in the lower panel 
of Figure is mass variation across sequences. The models 
shown in Figure [^confirm this. If we take the example given 
above for the ridge at (log P1[Wjk=12], logPi/Pl) « (1.3, 
0.3), panel A of Figure]^ shows that the period ratio of ra¬ 
dial order 2 to radial order 4 modes decreases as the model 
mass increases, just as suggested by the slope of the ridges. 

It was noted in Section that the ratio of the periods 
of the 1=1 and £ = 0 modes of radial order 4 increased 
slightly from the long period side to the short period side 


^ [Takayama, Saio, &: Ita| ( |2013^ use a nomenclature in which the 
radial order is offset by one from the nomenclature used here. For 
example, their p-mode P 2 corresponds to radial order 3 in this 
paper. 
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Figure 10. Plots of Wjk against the primary period PI for stars 
in the OGLE III catalogue (black points). The coloured lines show 
the positions of modes from the first five radial orders for stars 
of mass 0.8, 1.6 and 2.4 M0. The lines for each mode are plotted 
only in regions where they cross the sequence in which the mode 
is expected to dominate. The first two radial orders are shown 
as continuous lines and they correspond to the fundamental and 
first overtone radial pulsation modes. The lines for the next three 
radial orders are shown as dotted lines because the periods shown 
have been increased by AlogP = 0.043 from the period of the 
radial modes in an attempt to simulate the periods of the 1=1 
pulsation mode. The colours can be seen in the online version of 
the paper. 


of sequence A (the opposite slope to the cases considered in 
the previous paragraph). This result provides evidence for 
the mass dependence of the period ratio of the 1=1 and 
£ = 0 modes, with the period ratio increasing with mass. 
Gombining the results for the variation of the period ratio 
across sequence A from Section [3.2| with the results for the 
variation of mass across sequence A from Section |5.2| we 
find that dlog{P{£ = 1)/P{£ = 0))/dlogM « 0.036 for the 
stars on sequence A. This result provides a test for nonradial 
pulsation models of luminous red giant stars. 


7.3 Pulsation masses from the models 

Until now in this paper, only period ratios have been com¬ 
pared to models because period ratios are less sensitive to 
model parameters than are the periods themselves. In ad¬ 
dition, to compare model and observed periods in (say) the 
{Wjk, log P) plane, conversions from model Teff and L to the 
observed Wjk are required. These conversions are not par¬ 
ticularly well known for cool stars such as LPVs and there 
are additional problem due to interstellar and circumstellar 
reddening and the distance modulus (although the latter is 
now well-determined for the LMG). Keeping these problems 
in mind, in this final sub-section, we do compare the models 
to the observations in the (Wjk, logP) plane. 

For each stellar model, Teff and L have been converted 
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to J and K in the 2MASS system using the tabulated bolo- 
metric corrections obtained from the model atmosphere cal¬ 
culation by Houdashelt et al. 120001 and|Houdashelt, Bell, fc 


Sweigart (20001 for a metal abundance [Fe/H] = —0.5 which 


is appropriate for the LMC and which is consistent with the 
metal abundance used in the model calculations. Conversion 
of J and K in the system of Houdashelt et al. (20001 and 


Houdashelt, Bell, fc Sweigart 120001 to the 2MASS system 
is done using the conversions in Carpenter ( |2001[ ). An LMC 
distance modulus of 18.54 and a reddening E{B — V) of 
0.08 from Keller & Wood (20061 are used, along with the 
reddening law of Cardelli, Clayton, & Mathis (19891 to get 
the reddening at J and K. 

Figure 10 is a plot in the (Wjk, logP) plane of the 
periods and luminosities of modes from the first five radial 
orders of stars with masses of 0.8, 1.6 and 2.4 Mq. For each 
mode, only a short line segment is plotted centred on the 
position where it crosses the sequence in which the mode is 
assumed to be dominant. Note that the models computed 
here only produce periods for radial modes. However, it was 
shown in Sections |4] and that the £ = 1 mode is dominant 
on sequences A', A and B, and in Section]^ it was shown 
that these sequences correspond to radial orders 5, 4 and 
3, respectively. Furthermore, in Sections |3.2| and it was 
shown that the period of the £ = 1 mode is longer than 
the period of the £ = 0 mode by AlogP « 0.043. Hence, 
to compare the models with the observations, the computed 
periods for radial orders 3, 4 and 5 have been increased 
by AlogP = 0.043 before being plotted in FigureThis 
correction is made in order to simulate the £ = 1 mode 
periods. The resulting periods are shown in Figure by 
dashed lines. 

We now make a detailed comparison of the models with 
the observations. Looking at sequence C in Figure |10[ it 
can be seen that fundamental mode pulsation in stars with 
masses in the range 0.8-2.4 M© explains most of these stars, 
with some stars outside this mass range required to explain 
the full length of the sequence (note that the present-day 
population of red giant stars in the LMC peaks at masses 
around 1.6 Mq e.g. Nie, Wood, & Nicholls 20121. Stars in 
a similar mass range, but pulsating in the radial first over¬ 
tone mode can explain the stars on sequence CL Similarly, 
this mass range fits sequences B, A and A' but at lower lu¬ 
minosities, as expected from the empirical evolution scheme 
described in Section [5.11 The mass required to fit these three 
sequences at a given luminosity increases when moving from 
sequence B to A to A', consistent with the results derived 
in Section |5.2| from the empirical evolution scheme. In ad¬ 
dition, the width of sequence A requires a significant range 
of masses to be present at a given luminosity. This is also 
consistent with the results derived in Section |5j2l Given the 
various results presented in this paragraph, there appears to 
be reasonable agreement between the models and the obser¬ 
vations in the (VFjk, logP) plane in spite of the potential 
difficulties mentioned above. 

There is one parameter on which the models and the 
empirical evolution scheme are in relatively poor agreement. 
This is the transition luminosity at which stars move from 
one sequence to the next, especially in the case of the transi¬ 
tion from sequence C' to sequence C. Looking at the evolu¬ 
tionary track for a 1.6 M© star (i.e. the orange lines in Fig- 


101, we see that the first overtone mode does not cross 


the long period edge of sequence C' until Wjk ~ 9. Accord¬ 
ing to the empirical evolution scheme, the star should then 
jump to the short period edge of sequence C and start pul¬ 
sating in the fundamental mode at Wjk ~ 9. However, the 
period of the fundamental mode of the 1.6 M© star at this 
stage puts the star at the long period edge of sequence C or 
beyond. This suggests that the ratio of the fundamental to 
first overtone periods in the models is too large. Actually, 
the disagreement between the models and the observations 
is not as bad as just suggested. If we look at panel C' of Fig¬ 
ure 1^ we see that the logarithmic period difference A log P 
between the fundamental mode to the first overtone mode 
is ~0.26 around Wjk ~ 9 — 10. This is significantly larger 
than the difference AlogP ~ 0.18 between the long period 
edge of sequence C' and the short period edge of sequence C. 
At Wjk = 10, AlogP for the 1.6 M© star is indeed ~ 0.26 
but by Wjk = 9 AlogP is much too large at ~ 0.41 so 
that there is still a discrepancy between the model period 
ratios for the first two modes and the observed ratio of these 
two modes. It is worth noting that the discrepancy discussed 
here seems to be worst at high masses and luminosities. At 
M — 0.8 Mq, the models predict a transition from the long 
period edge of sequence C' to the short period edge of se¬ 
quence C which is as predicted by the empirical evolution 
scheme. 

Finally, we take a look at the period difference between 
RGB and AGB stars. The lines drawn in Figurej^belong to 
both RGB and AGB stars. The computed separation of the 
RGB and AGB periods for each mode at a given luminosity 
is almost imperceptible on the plot (the separation can just 
be seen for the 0.8 Mq models). It is certainly much less 
than the period shift of AlogP ~ 0.05 evident on sequence 
A when moving from the RGB at Wjk > 11.3 to the AGB at 
Wjk < 11.3. The observed period shift can not be explained 
by mass loss between the RGB tip and the AGB since mass 
loss increases the period of pulsation whereas the period is 
observed to decrease. The models thus predict that there is 
a very small period change between the RGB to the AGB for 
a given mode. A possible interpretation of the much larger 
period shift observed on sequence A when going from the 
RGB to the AGB is that the AGB is dominated by radial 
(£ = 0) pulsators while the RGB is dominated by nonra- 
dial £ = 1 pulsators (remember that AlogP ~ 0.043 for the 
period ratio of the £ = 1 and £ = 0 modes). However, this 
result conflicts with the results in Section which suggested 
that the £ = 1 mode becomes more dominant at higher lu¬ 
minosity. The observational resolution of this conflict awaits 
larger data sets such as those that will come from OGLE IV. 


8 SUMMARY AND CONCLUSIONS 

Multiple pulsation modes in stars in the OGLE III catalogue 
of LPVs in the LMC have been examined. It has been shown 
that the ratios of the multiple periods in individual stars do 
not agree with the ratios of the periods of the sequences of 
LPVs in the PL diagram. This leads to the conclusion that 
the average masses of stars on the different sequences vary, 
with the shorter period sequences containing more massive 
stars. An examation of the evolution of stars in the PL di¬ 
agram shows that the mass variation between sequences is 
the same as the mass variation across individual sequences. 
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As an example, it is shown that the average mass of a star on 
sequence A is ~26% higher than the average mass of a star 
on sequence B and that the mass variation across sequence 
A is of a similar magnitude, with the more massive stars 
on the shorter period side of the sequence. In general, the 
masses of stars at a fixed luminosity on adjacent sequences 
are found to increase by ~ 16-26% towards the shorter pe¬ 
riod sequence. Fitting pulsation models to the period ratios 
in multimode pulsators shows that the sequences C, C', B, 
A and A' correspond to pulsation in radial orders 1 to 5, 
respectively (with radial order 1 corresponding to the fun¬ 
damental mode). Closely-spaced modes revealed by the pe¬ 
riod ratios in individual stars show the presence of nonradial 
modes. The observed period ratios indicate the presence of 
the ^ = 0, 1 and 2 modes on sequence A (radial order 4) and 
the i = 0 and 1 modes on sequence B (radial order 3). The 
relative populations and amplitudes of the modes of different 
I show that the 1=1 mode is the most common and of the 
largest amplitude in both sequences. The positions of sec¬ 
ondary mode sequences relative to primary mode sequences 
in the PL diagram also suggest that the 1=1 mode is the 
most prominent primary and secondary mode on sequences 
A', A and B and that the I = Q mode is the most prominent 
primary and secondary mode on sequences C' and C. This 
result is supported by the good ht of the period ratios in 
radial pulsation models to the observed period ratios. 
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